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Abstract 

Let S be a finite set of points in the Euciidean piane. Let D be a 
Delaunay triangulation of S. The stretch factor (also known as dilation 
or spanning ratio) of D is the maximum ratio, among all points p and q 
in S, of the shortest path distance from p to q in D over the Euclidean 
distance | \pq\\. Proving a tight bound on the stretch factor of the Delaunay 
triangulation has been a long standing open problem in computational 
geometry. 

In this paper we prove that the stretch factor of the Delaunay triangu- 
lation of a set of points in the plane is less than p = 1.998, improving the 
previous best upper bound of 2.42 by Keil and Gutwin (1989). Our bound 
1.998 is better than the current upper bound of 2.33 for the special case 
when the point set is in convex position by Cui, Kanj and Xia (2009). This 
upper bound breaks the barrier 2, which is significant because previously 
no family of plane graphs was known to have a stretch factor guaranteed 
to be less than 2 on any set of points. 



*A preliminary version of the paper appears in the Proceedings of the 27th Annual Sym- 
posium on Computational Geometry (SoCG 2011). 
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1 Introduction 



Let S be a finite set of points in the Euclidean plane. A Delaunay triangulation 
of S is a triangulation in which the circumscribed circle of every triangle contains 
no point of S in its interior. An alternative equivalent definition is: An edge xy 
is in the Delaunay triangulation of S if and only if there exists a circle through 
points x and y whose interior is devoid of points of S. A Delaunay triangulation 
of S is the dual graph of the Voronoi diagram of S. 

Let D be a Delaunay triangulation of S. The stretch factor (also known as 
dilation or spanning ratio) of D is the maximum ratio, among all points p and q 
in S, of the shortest path distance from p to q in D over the Euclidean distance 

I Im- 
proving a tight bound on the stretch factor of the Delaunay triangulation 
has been a long standing open problem in computational geometry. Chew [5] 
showed a lower bound of tt/2 on the stretch factor of the Delaunay triangulation. 
This lower bound of tt/2 was widely believed to be tight until recently (2009) 
when Bose et al. [2] proved that the lower bound is at least 1.5846 > 7r/2, which 
was further improved to 1.5907 by Xia and Zhang [T7]. Dobkin, Friedman, 
and Supowit [7J [5] in 1987 showed that the Delaunay triangulation has stretch 
factor at most (1 + \/5)n/2 « 5.08. This upper bound was improved by Keil 
and Gutwin [H [13] in 1989 to 2?r/(3 cos (vr/6)) w 2.42, which currently stands 
as the best upper bound on the stretch factor of the Delaunay triangulation. 
For the special case when the point set S is in convex position, Cui, Kanj and 
Xia [6] recently proved that the Delaunay triangulation of S has stretch factor 
at most 2.33. 

In addition to its theoretical interest, improving the upper bound on the 
stretch factor of the Delaunay triangulation has a direct impact on the problem 
of constructing geometric spanners of Euclidean graphs, which has significant 
applications in the area of wireless computing (for more details, see |15j). Many 
spanner constructions in the literature rely on extracting subgraphs of the De- 
launay triangulation (see for example [3) [9j [TOj [14]) and their spanning ratio 
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is expressed as a function of the stretch factor of the Delaunay triangulation. 
Hence improving the upper bound on the stretch factor of the Delaunay trian- 
gulation automatically improves the upper bounds on the spanning ratio of all 
such spanners. 

In this paper we show that the stretch factor of the Delaunay triangulation 
of a point set in the plane is less than p — 1.998, improving the current best 
upper bound of 2.42 by Keil and Gutwin jT2] H3] • Our bound 1.998 is better 
than the current upper bound of 2.33 for the special case when the point set is 
in convex position [6]. 

The notion of stretch factor can be defined on any plane graph G as the 
maximum ratio, among all vertices u and v in G, of the shortest path distance 
between u and v in G over the Euclidean distance The family of plane 

graphs with the best known upper bound on the stretch factor on any set of 
points has a stretch factor of 2 [H[S]. Our result improves this bound to 1.998. 

Our approach in proving the upper bound on the stretch factor of the De- 
launay triangulation is significantly different from the previous approaches 
[H [TSJ, [T3l [6] . Our approach focuses on the geometry of a "chain" of circles in 
the plane. The "stretch factor" of a chain can be defined in analogy to that of 
the Delaunay triangulation. We prove an upper bound on the stretch factor of 
a chain and derive the same upper bound on the stretch factor of the Delaunay 
triangulation as a special case. 

The paper is organized as follows. The necessary definitions are given in 
Section [2] The main theorem of the paper is given in Section [3] There we 
also discuss the proof strategy and provide an outline of the proof for the main 
theorem. Most of the technical details are captured by two lemmas, whose 
proofs appear in Section 2] and Section [SJ The paper ends with a discussion on 
the possible improvements of our approach in Section [SJ 
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2 Preliminaries 



We label the points in the plane by lower case letters, such as p, q, u, v, etc. For 
any two points p, q in the plane, denote by pq a line in the plane passing through 
p and q, by pq the line segment connecting p and q, and by pt[ the ray from p to 
q. The Euclidean distance between p and q is denoted by \\pq\\. The length of a 
path P in the plane is denoted by \P\. Any angle denoted by Zpoq is measured 
from to in the counterclockwise direction. Unless otherwise specified, the 
angles are defined in the range (— it, n]. 

Definition 1. We say that a sequence of distinct finite circlefQo = (Oi, O2, ■ ■ ■ , O n ) 
in the plane is a chain if it has the following two properties. Property (1): 
Every two consecutive circles Oj,Oj+i intersect, 1 < i < n — 1. Let dj and &i 
be the shared points on their boundary (in the special case where Oi,Oi+i are 
tangent, a% = hi). Without loss of generality, assume aj's are on one side of 

and bi's are on the other side. Denote by c[ l the arc on the boundary 

(i) 

of Oi that is in Oi+i, and by C7Vj the arc on the boundary of Oi+i that is in 
Oi. We refer to c|' +1 ^ and C^_\ as "connecting arcs". Property (2): The 
connecting arcs cj 1- ^ and C t - ^ on the same circle Oi do not overlap, for 
2 < i < n — 1; i.e., C^ 1 ^ and Cj* +1 ' do not share any point other than a 
boundary point. See Figure [T] for an illustration. Denote by Oij a sub-chain of 
O: Oij = '() (),>. 

Given a chain O = (0±, O2, ■ ■ ■ , O n ). Let u be a point on the boundary of 
0\ that is not in the interior of O2 ■ Let v be a point on the boundary of O n that 
is not in the interior of O n -\. We call u,v a pair of terminal points (or simply 
terminals) of the chain O. For notational convenience, define oq = bo = u and 
a„ = b n = v. Every circle Oi has two arcs on its boundary between the line 
segments ai-ibi-i and a;&;, denoted by Ai and Bi. Without loss of generality, 
assume that ai_i,Oj are the ends of Ai and are the ends of Bi, for 

1 < i < n. This means that Ai . . . A n is a path from « to d on one side of the 

1 In this paper, a circle is considered to be a closed disk in the plane 
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Figure 1: A chain O. The connecting arcs are green (gray in black and white 
printing). The connecting arcs on the boundary of the same circle are disjoint. 
In this example, the terminals u, v are obstructed while the terminals u, v' are 
unobstructed. 

chain and B\ . . . B n is a path from u to v on the other side of the chain. An 
arc Ai or Bi may degenerate to a point, in which case a.;_i = a, or 6j_i = bi, 
respectively. 

Let Dq(u, v) = upi . . .p n ^iv be the shortest polyline from wtow that con- 
sists of line segments upT, PiP2, ■ ■ ■ ,Pn-iv, where pi € a^bi for 1 < i < n— 1. The 
length of Da(u,v), denoted by \Do(u,v)\, is the sum of the lengths of the line 
segments in Dq{u, v). See FigureQ] If the polyline Do{u, v) contains a point pj 
which is a,j or bj, for some 1 < j < n— 1, we say that u, v are obstructed. If u, v 
are unobstructed, then Dq{u, v) is the line segment uv. We have the following 
fact. 

Fact 1. If Do(u,v) is the line segment uv, then ui stabs 0\, . . . , O n in order. 
That is, for any 1 < i < j < n, Uv enters Oi no later than entering Oj and 
exits Oi no later than exiting Oj . 

Proof. For 2 < i < n — 1, the two connecting arcs on Oi do not overlap. 
Hence a,;_i&i_i and ai + \bi + \ must appear on the different sides of the line a,6j. 
Without loss of generality, assume that aibi is a vertical line which divides the 
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plane into two half-planes, ai_i&;_i appears in the left half-plane and aj+ifej+i 
appears in the right half-plane. If Do(u,v) is the line segment uv, then uv 
intersects from left to right. Thus ui enters Oj_ i and Oi in the left half- 
plane. Since Oi is contained in Oi_i in the left half-plane, uv must enter Oi—± 
no later than entering Oi. Likewise, uv exits Oi-i and Oi in the right half-plane. 
Since Oi-\ is contained in Oi in the right half-plane, uv must exit Oi-i no later 
than exiting Oi. This completes the proof. □ 

The polyline Dq{u,v) can be visualized as a "rubber band" connecting u 
and v "inside" the circles in O. The polyline Dq{u,v) can be regarded as 
a correspondence of the straight line between two point u, v in the Delaunay 
triangulation. Note that for the purpose of bounding the stretch factor of the 
Delaunay triangulation, only the case where u, v are unobstructed is relevant. 
However, for our proof to work, it is necessary to consider the case when u, v 
are obstructed (see Section [4} . 

We define the shortest path between u and v in O, denoted by Po(u,v), to 
be the shortest path from u to v while traveling along arcs in {Ai, . . . ,A n } U 
{B±, . . . , B n } and line segments in {a\b\, . . . , a n _i& n _x}. Note that the vertex 
set of Pq{u, v) is a subset of {ao, . . . , a„} U {bo, . . . , b n } and the edge set of 
Po{u,v) is a subset of {Ai, . . . ,A n } U {B x , .. ., B n } U {ai&i, . . . , a„_i6„_i}. Its 
length, \Pq{u, v)\, is the total length of the edges in Pq{u,v). For example, in 
Figure [U Pq(u,v) is the shortest path from u to v while traveling along the 
thick arcs and lines. Pq(u,v) corresponds to the shortest path between two 
point u,v in the Delaunay triangulation. 

Now we can define the stretch factor of a chain O, denoted by Co, as the 
maximum value of 

\Po(u,v)\/\D (u,v)\, (1) 

over all terminals u, v of O. The stretch factor of a chain is analogous to that 
of a Delaunay triangulation. 

In this paper, we will focus our attention on the stretch factor of a chain. 
Compared to the Delaunay triangulation, it is much easier to manipulate a 
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chain. For example, it is easy to transform a single circle in a chain without 
affecting other circles in the chain. In contrast, changing a triangle in a Delaunay 
triangulation always affects other triangles. 

3 The Main Theorem 

The following is the main theorem of this paper, which gives an upper bound 
on the stretch factor of a chain. 

Theorem 1. For all O, the stretch factor Co is less than p, where p — 1.998. 

Assuming Theorem [JJ is true, we can derive an improved upper bound on 
the stretch factor of the Delaunay triangulation as a special case. 

Corollary 1. The stretch factor of a Delaunay triangulation D of a set of points 
S in the plane is less than p — 1.998. 

Proof. For any two points x,y € S. Let T be the sequence of triangles in D 
crossed by x%. Let O be the corresponding sequence of circumscribed circles 
of the triangles in T. It is clear that O is a chain and x, y are terminals of 
O. The shortest path distance from x to y is at most \Po(x,y)\ and \\xy\\ — 
\Do(x,y)\. By Theorem [TJ the stretch factor of a Delaunay triangulation is at 
most \P (x,y)\/\\xy\\ = \P {x,y)\/\D (x,y)\ < p. □ 

Due to the complex nature of the proof for Theorem [TJ in this section we 
discuss the proof strategy and present a outline of the proof. The main technical 
details of the proof are captured in two lemmas, whose proofs are given in the 
subsequent sections. 

When O has only one circle, it is easy to see that for all O, it, and v, 
\Po(u, v)\/\Dq(u, v)\ < 7r/2 < p. So it is natural to attempt an inductive 
proof based on the number of circles in O. A simple induction would require 
us to show that adding a circle to a chain will not increase the stretch factor. 
However, this does not work because one can always increase the stretch factor 
of a chain by adding a circle to it, albeit by a very small amount [17]. We 
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tackle this problem by amortized analysis. Specifically, we introduce a potential 
function $e> and define a target function To{u,v): 

T (u,v) = \Pq(u,v)\ - X\D (u,v)\ +$0, (2) 

where A = 1.8 is a parameter whose value is determined by the potential func- 
tion. The goal is to prove Theorem[T]by showing To{u,v) < for all 0,u, and 
v. 

The key component of our strategy is the selection of an appropriate poten- 
tial function <£>e>, which is described in the following. 

Definition 2. Let Oj_i and Oi be two consecutive circles in a chain O. Without 
loss of generality, assume that their centers Oj_i,Oj lie on a horizontal line and 
that CLi-i is on or above the line Oi-\Oi. Refer to Figure [21 Let q~L^ be the 
"peak" of Oi—i with regard to Oi-iOi, i.e., the point on the upper boundary of 
0,_i that is farthest from the line Oj-io,-. Likewise, let q*~ be the "peak" of Oi 
with regard to Oi-\Oi (the sign or *~ indicates whether the peak is defined 
with the previous circle or the next circle in O). Let Q~Li be the upper arc 
between q^_ l and a^_i on the boundary of Oj_j and let be the upper arc 
between q^~ and a%-i on the boundary of O;. If Q~Xi is inside Oi, we color it red; 
otherwise we color Q~L 1 green. Likewise, we color Q*f red or green depending 
on whether it is inside Oi-\. Let Vi be a path from q~li to q£~ consisting of 
Q^Li and Q'f . Let Hi be the horizontal distance traveled along the path Vi with 
green arcs contributing positively to Hi and red arcs contributing negatively to 
Hi. Similarly, let Vi be the vertical distance traveled along the path Vi with 
green arcs contributing positively to Vi and red arcs contributing negatively to 
Vi. The potential function is defined as follows: 

n 

$o = <p(r n -r 1 )-^y2(2H i + V i ), (3) 

where r» is the radius of Oi and ip = -^(1 — X/p) is a parameter that determines 
the "weight" of the potential function. 

We have the following fact. 



8 




Figure 2: The definition of Vi — the thick colored arcs (in black and white print- 
ing, red is the darker color and green is the lighter color). In (a) and (b) V% 
has one green arc and one red arc. In (c) Vi has two green arcs, (c) shows an 
example of V% and Vi+i defined on three consecutive circles Oi-\, Oi and Oi+i. 
Note that Oi has two different peaks. The peak q~ is defined with the previous 
circle Oj_i. The peak is defined with the next circle Oi+i. 
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Fact 2. Let r.i and rj_i be the radii of Oi and Oi-\, respectively. Then Hi = 
\\oiOi-\\\ and Vi > \n - r*_i|. 

Proof. For this proof, fix a coordinate system where the origin is Oj, a;-axis is 
Oi-io\, and a%-\ is above the x-axis. Let X q r+_ and y g r» be the x- and y- 
coordinates of q~l 1 . Let -Xq<- and Y q *- be the a;- and y-coordinates of q£~. Let 
X ail and Y «_i be the x- and y-coordinates of Ot_ i. We distinguish three cases. 

• Case 1. QX'-\ is green and is red. See Figure Ufa) . Then 

Hi = (X ai l — Xqr*^ — (X ai l — X q *-) — X q < X q r^ i — 1 1 0« O-j — 1 j ] . 

In this case, < rj_i. We have 

Vi = ( Y q7li - Y ai _ t ) - (Y qt - Y ai _ t ) > Y q ^ - Y qt = r;_i - r< = \n - r 4 _i| . 

• Case 2. is red and Q*~ is green. See Figure [2fb) . Then 

-Hi = -^Oj-i) — {Xq->_ t — ^ ai _i) = -^gj ^9^1 = 1 1 °i °i — 1 1 1 - 

In this case, n > rj_i. We have 

Vi = {Y qt -Y ai _,)-{Y q r ti -Y ai _ 1 )>Y qt -Y q r ti =r i -r i - 1 = |rj-ri_i|. 

• Case 3. Both Q^_ 1 and are green. See Figure (2fc). Then 

Hi = (X qt - X ai _J + {X ai _, - X^J = X qt - X q ^ = \\o i0i -x% and 

V = (Y qr - r aj _j + (y^ - y^j > |y ?r - y^j = \ n - r <_i|. 

So in all cases, the statement is true. This completes the proof. □ 

The potential function is designed with three goals in mind. 
First, the potential function $e> is designed such that adding a circle to O 
does not increase $e>, as shown below. 

Proposition 1. $e> < $Oi,„-i- 
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Proof. By the triangle inequality, ||o n o„_i|| > ||o n a„-i|| - ||o n _ia n _i|| = r„ - 
r n -\. Combining this with Fact[2J we have 

2H n + V n > 211 OnOn — 1|| ~r \r n ^n— 1 

> 2(r„ - r„_i) + (r„ - r n _i) 

= 3(r» - r„_i). (4) 

So $ - = ^(r„ - r„_i) - f (2tf„ + 14) < 0. □ 

Secondly, the potential function $q is designed such that Ta{u,v) < for 
all O, u, and v. This is Lemma [TJ the main technical lemma of the paper, whose 
proof is given in Section [4] 

Lemma 1. For all O, u, and v, Tq(u, v) < 0. 

Thirdly, the potential function $o is designed such that its value can be 
bounded from below as a function of |Pe>(u,u)| for some chain O with certain 
extremal properties. This is Lemma [H whose proof appears in Section [5] 

Lemma 2. Let O be a set of chains whose stretch factor is greater than or equal 
to a threshold t. If O is not empty, then there exists a chain O* € O with ter- 
minals it, v such that \Po* (u, v)\/\Dq* (it, v)\ > r and <&q* > ~^^-\Pq' (u-, v)\. 

With Lemma Q] and Lemma [5J we are ready to give the proof for the main 
theorem. 

(proof of Theorem^. For sake of contradiction, suppose that there is a non- 
empty set O of chains O with stretch factor Co > P- By Lemma [2] there exists 
a chain O* € O with terminals u and i> such that 

\Po*(u,v)\/\D *(u,v)\>p, (5) 

and 

$ .>_^|P .( ti)1 ,)|. ( 6 ) 
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By Lemma [TJ 

To.(«,tj) = \P *(u,v)\-\\D *(u,v)\ + $o* <0. (7) 
Combining ^ and ([7]), we have 

\Po-(u,v)\ - \\Do*(u,v)\ - ^.\p .(u,v)\ 
< |P . (u,t;)| - \\D *(u,v)\ + $o* 

= T *(u,v)<Q. (8) 
Recall that <p = -^=(1 — X/p). Rearranging ©, we have 

\Po*(u,v)\ < ^\D *{u,v)\ 



1 3 



-\D .(u,v)\ 



1 - (1 - A/p) 1 
= p\Da*(u,v)\. (9) 

This is a contradiction to ([5]). Therefore Co < p for all O. □ 

The rest of the paper contains the proofs for Lemma [1] and Lemma [2] 

4 Proof of Lemma 1 

This section is devoted to the proof of Lemma [TJ 

(proof of Lemma{^l. Proceed by induction on n, the number of circles in O. 
If ri = 1, then <£>£> = 0, Po(u,v) is the shorter arc between u and v on the 
boundary of 0\, and Do{u,v) is the chord between u and v inside 0\. So 
\Po(u,v)\ < n/2\Do(u,v)\ and hence To(u,v) = \Po(u,v)\ — X\Do(u,v)\ + 
$o = \Po(u,v)\-X\D (u,v)\ < 0. 

Now assuming that the statement is true when there are less than n circles 
in O, we will prove that it is true when there are n circles in O. 

For the rest of this section, fix an arbitrary chain O = (Oi, O2, ■ ■ ■ , O n ) and 
an arbitrary terminal u on the boundary of 0\ . We start by considering some 
special cases. 
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First consider the case when v is either a„_i or 6 n _i. 
Proposition 2. To(u, a n -i) < and Te>(u, 6 n _i) < 0. 

Proof. Consider the sub-chain 0\ ;n -i = (0±, . . . , 0„_i). Since u, a n _i are ter- 
minals of Oi^ n -i, by the inductive hypothesis, 

f Oi,„_i(«)On-l) 

= l-PGi.n-i ( w i a n-l)l - A|Doi,„_i («, On-l)| + *Oi,„-i 

< 0. (10) 

It is easy to verify that \Pq(u, a„-i)| < l-Poi „-i( u > a n-i)| an d |Z?e>( u 5 a n-i)| = 
\Do 1 n _ 1 (u, a n —i)\. Combining these with (|10[) . we have 

To(u,a„_i) = \Po(u,a n -i)\ - X\Do(u, a n _i)| + $e> 

< |-Po,, n -i( u > a n-i)| - A|£> 0l i „_ 1 (u, a„_i) | + $o 

< T c , l 7l _ 1 (u,a„_i) + $o - 

<0. (11) 

The last inequality is from Proposition [T] Similarly, we have To(u, &„-i) < 
0. □ 

Next consider the case when u and v are obstructed. 
Proposition 3. If u and v are obstructed, then To(u,v) < 0. 

Proof. If u and v are obstructed, then Dq(u,v) contains a point pj that is 
either <ij or bj, for some 1 < j < n — 1. Without loss of generality, as- 
sume pj — aj. Consider two sub-chains of O: Oi.j+i = (Oi, . . . , Oj+i) and 
&j+i,n — {Oj+i, ■ • • , O n ). Points u, aj are terminals of Oij+i and points a-,-, v 
are terminals of Oj+i j7l . If j = n — 1, by Proposition [5J 

Toi, j+ i(u,aj) = T («,o n _i) < 0. (12) 
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If j < n — 1, has less than n circles, and by the inductive hypothesis 

To 1J+1 (u,a 3 ) < 0. (13) 

On the other hand, Oj+i >n has less than n circles, and by the inductive hypoth- 
esis, 

To J+1 , B (ai,«)<0. (14) 

Note that 

\Po(u,v)\ < \P 0l , ]+1 {u,a 3 )\ + \P 0]+1 .M,v)l (15) 
\D (u,v)\ = \Do 1 , j+1 (u,a j )\ + \D 0j+1 Ja j) v)\, (16) 

and 

$0 = $e>i a+1 +$o i+1 ,„. (17) 
Combining (fT2|) — (fT7|) . we have 

T (u,») = |Po(«,«)| - A|D («,«)| + $c 

< l-Fbi, 3+ i - A l- D d, 3+ i I + *Oi, J+ i 

+ |Poi+i.»(°J> w )l - x \Do J+1 ,Ja 3 ,v)\+<5>o 3+1 , n 
= ^o ltj+1 (u,aj) + T e . 3 . +1 _„(a i ,w) 

< 0, (18) 
as desired. □ 

Now we only need to consider the set of unobstructed terminal points v on 
O n , which form an arc, denoted by A and referred to as the unobstructed arj^j. 
Denote by P@ " (it, u) and Pq ™ (u, u) the shortest paths from u to v that includes 
A„ and B n (on O n ), respectively. We call v pivotal if \Pq" (u, v)\ = \P@ n (u, v)\. 



2 For notational convenience, assume that A includes its boundary points. This makes A a 
closed set. 
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Figure 3: Illustration for Proposition HI The angle Zuvo n is currently in range 



below uo„, Zuvon becomes negative. 

Proposition 4. The maximum of Tq(u,v) occurs when v is on the boundary 
of A or when v is pivotal in A. 

Proof. Let A' be an arbitrary sub-arc of A that does not contain a pivotal point 
in its interior. In order to prove the proposition, it suffices to show that the 
maximum of Tq(u, v) for all v G A' occurs when v is on the boundary of A' . 

Since A' does not contain a pivotal point in its interior, either \Pq" (u, v)\ < 
I-Pq™ {u,v)\ for every point v E A' or \P^ n (u,v)\ < \Pq" (u,v)\ for every point 
v e A'. Without loss of generality assume that \Pg n (u,v)\ < \Pg n (u,v)\ for 
every point v in A'. Fixing other parameters, To(u,v) is a function of \A n \ 
as v moves along A' . One observes the following. First, remains constant 
when v moves along A' . Secondly, \Pq{u, v)\ is a linear function of \A n \ because 
\Po(u,v)\ = \P n (u,v)\ = \Po(u,a n -i)\ + \An\, where \Po(u, a n -i)| remains 
constant when v moves along A' . Thirdly, — \\Dq{u, v)\ is a convex function of 
\A n \. To see why, refer to Figure [3] One observes that 



(0, it/2). As v moves away from <z n _i 



along A', Zuvo n decreases. When v moves 



d\D (u,v)\ 
d\A n \ 



cos a = sin(Zuvo n ). 



(19) 
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Also observe that Zuvo n decreases as v moves away from a n -i along A' and 
hence < 0. Since v is the exit point of ray m) on O n , Zuvo n G 

(-tt/2, tt/2) and hence dsi f/ uv °^ > . So dsil \( z ""°") = ^in( Zu TO „) < 

v fill dZuvo n — a|A^| dZuvo n d\A n \ — 

0. Combining this with (|19p . we have 

d 2 (-X\D ju,v)\) _ d(-Xsm(Zuvo n )) 

d\A n \i ~ d\A n \ ~ ' 

which implies that — X\Do(u, v)\ is a convex function of \A n \. 

Being a sum of convex functions, Tq(m, = \Po(u, v)\ — X\Dq(u, v)\ + &o 
is also a convex function of \A n \. This proves that the maximum of To(u,v) 
appears on the boundary of A', as desired. □ 

Therefore, if we can show that To(u, v) < when v is on the boundary of 
the A and when v is pivotal, then To (it, v) < for any point v in A and we are 
done. 

Proposition 5. Tq{u,v) < if v is on the boundary of A. 

Proof. If v — a n -i or v = b n -i, by Proposition [21 To(u,v) < 0. If v is 
neither a n _i nor b n -i, then u and u must be obstructed and by Proposition [3l 
T o (u,v)<0. □ 

What remains to be shown is that Tq(u,v) < when v is pivotal. This 
requires a careful analysis of the geometry. To start, we fix a coordinate system 



where the origin is the center point of a n -ib n -i, the x-axis is o„_io n where 
o n —i is to the left of o„, and the y-axis is a„_i6„_i where a„_i is above b n -\- 
See Figure |4]for an illustration. Since uv crosses a n -ib n -i from left to right, u 
is to the left of the j/-axis and v is to the right of the y-axis. By flipping the 
geometry along the a:-axis, if necessary, we can assume Y v < Y u , where Y v and 
Y u are the y-coordinates of v and u, respectively. Let q be the exit point of the 
x-axis on O n . Let v' be the entry-point of uti on O n . 

We define the following parameters. Let a = Zqo n a n -\ and ft = Zvo n q. Let 
7 be the angle from ui to the rr-axis in the counterclockwise direction. Their 
ranges are given as follows. 
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Figure 4: Illustration of the coordinate system, the definition of the parameters 
a, P, 7, and the movement of o n toward o n —\. 

Since a n _i is on or above the x-axis, < a < iz. If a = or a = it, then 
a n _i = 6„_i, which means that D&(u,v) contains a„_i and Tq(u,v) < by 
Proposition [3] So we can assume 

< a < ir. (20) 

Since v is a pivotal point, \Pq" (u,v)\ = \P@ n (u,v)\, where \P^ n (u,v)\ = 
|Po(«,an_i)| + \A n \ and = 

\Po(u,b n -i)\ + \B n \. Since a n _i and 6„_i are connected by a line segment 
a n _i6 n _i, wehave-||a n _i& n _i|| < \Po(u,a n -i)\-\Po(u,b n -i)\ < ||a„_i6„_i||. 
So -||a n _i6„_i|| < \A n \ - \B n \ < ||a„_i&„_i||, where ||a„_i6 n _i|| = 2r„sina 
and \A n \ — \B n \ — 2r n f3. Thus the range of j3 is 

— sin a < (3 < sin a. (21) 

We have < 7 because Y v < Y u . Since uv crosses a„_i6 n _i, the largest value 
of 7 occurs when uv passes through a n — \. This means 7 < 7r/2 — Z.b n —\a n —xv — 
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7r/2 — (a — (3)/2, since Z&„_ia Tl _iv = (a — /3)/2. So the range of 7 is 

< 7 < tt/2 - (a - /3)/2 < tt/2. (22) 

The last inequality is true because a — /3 > a — sin a > 0. 

We proceed by distinguishing two cases depending on the value of 7 in com- 
parison to a threshold value 7 + defined as 

, 3 sin a — a ( a + sin a \ , . 

7 = : 1- arcsin rrr- — r • (23 

Proposition 6. If v is a pivotal point and 7 > 7 + , then To(«,«) < 0. 

Proof. Since u, v are unobstructed, ui crosses all line segments ajfej, 1 < i < 
n — 1. The points {ai, . . . , a„_i} are all above the line uv. Recall that Do(u, v) 
is a "rubber band" between u and u. When one end of the "rubber band" moves 
from v to a„_i along the boundary of O n , the result is polyline Do(u,a n -i), 
which is a path from u to a„_i that is convex-away from ua n -\ and is between 
ua„_i and ui. 

See Figure 131 Let a/j be the last turning point in the polyline Do(u,a n -i). 
Then the part of Dq{u 1 a n _i) between and a„_i is a line segment, and 
afc, a„_i are terminals of the sub-chain Ok+i, n - By Fact[TJ afca„_i stabs Ok+i , ■ • • , O n 
in order. Since 7^ afe + i (otherwise dfc+i would be the last turning point), 
is an entry point on Ok+i- This means that appears no later than the entry 
point of Do(u,a n -\) on O n . Combining this with the fact that Dq(u, a„_i) 
is between ua„-\ and uv, we know that at- is in the triangle Aw'a n _i. Now 
ua n -i U Dq{u, a n -i) is a convex polygon which contains afca n _i as an edge. 
So the whole path Do(u,a n -i) is above the line ossein— 1. This implies that 
Dq(u, a n —i) lies in the triangle AWa n _i. See FigureU]for an illustration. 

Since Do(u, a n _i) is a path from it to a n _i that is convex-away from ua„_i 
in the triangle A.uv'a n —i- By convexity, we have |Do(it, a„_i)| < + 
||u'a„_i|| (see [TJ p. 42]). Recall that \Dq(u,v)\ = So 

\D (u,v)\ - \D (u,a n -i)\ > \\uv\\ - (\\uv'\\ + ||u'a n _i||) 

= \\v'v\\ — \\v'a n -x\\. (24) 
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Refer to Figure^ \\v'v\\ — 2r n cos(Zvv'o n ) and ||i/a n _i|| = 2r n cos(Zo n v'a n -i), 
where Z.vv'o n = Zo n vv' = — 7 and Zo n v'a n —i — Zvv'a n —i — Zvv'o n — 
(a + 13)/ 2 - (0 - 7) = a/2 - 0/2 + 7. From (gl]) and by the trigonometric 
identity, we have 

\Do(u,v)\ - \D (u,a n -i)\ 

> \\v'v\\ - \\v'a n -l\\ 

= 2r„ cos(0 - 7) - 2r„ cos(a/2 - (3/2 + 7) 

= -4r„ sin(a/4 + 0/4) sin(30/4 - a/4 - 7). (25) 

Since v is pivotal, \Pq{u, v)\ = \Po{u, a„_i)| + \A n \, where \A n \ = r n (a + (3). 
By Proposition [5J Te>(u,a„_i) < 0. Combining these with (|23|) . we have 

T (u,v) = \P (u,v)\-\\D {u,v)\+$o 

= |P (u,a»_i)| + \A n \- X\D (u,v)\ + $ 

= T («,o„_i) + \A n \ - \(\D (u,v)\ - |Do(u,o n _i)|) 

< \A n \ - X(\D (u,v)\ - |D (u,a n -i)|) 

< r„(a + 0) + 4r„Asin(a/4 + ,3/4) sin(30/4 - a/4 - 7). (26) 
Define a function 

h{a, 13, 7) = r„(a + 0) + 4r n A sin(a/4 + 0/4) sin(30/4 - a/4 - 7). 

Then 

— = -4r„Asin(a/4 + 0/4)cos(30/4-a/4-7). 
07 

By Ipljl. we have a/4+0/4 > a/4-sina/4 > 0, and a/4+0/4 < a/4+sina/4 < 
tt/4. So 

< a/4 + 0/4 < tt/4. (27) 

By {25]), we have 30/4 -a/4 -7 > 30/4- a/4 - tt/2 + (a - 0)/2 = a/4 + 0/4 - 
tt/2 > -tt/2. Also 30/4 - a/4 - 7 < 30/4 < 3/4. So 

-tt/2 < 30/4- a/4-7 < 3/4. (28) 
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From ||27H and fl2SJ), sin(a/4 + /3/4) > and cos(3/?/4 - a/4 - 7) > 0, and hence 



Let 

7*(a,fl) = 33/4- a/4 + arcsin ( 1 • 

' V M/ 7 \4Asin(a/4 + /3/4)/ 

Then h{a,(3,j*) = 0. Let /z = a/4 + ^/4 and f = By d2TJ) , < /z < tt/4. 

We have 



»M_*(a/4 + /*/4) >0> (3Q) 



9/3 9/3 
Since /it < tan// and sin/z > for < /z < tt/4, we have 

9f _ ^XsinTI 



(1 - /V tan /z)/(A sin /ti) > 0. (31) 



dfi dpi 

From (Ell), i/ = ir ^— < , Z^ 4 .,, < 0.618. Also v > 0. So 

1 !" Asin/i — Asin(7r/4) 

9 arcsin v 
dv 

By idU), {31]), and (pi). 



= l/\/l-^ 2 >0, (32) 



97* . 9 arcsin^ 9^ 9/i 

95 = / + 9^ 9^' 9/3 > ' 

This means 7* is increasing in /3, and hence 

7* (a, /3) < 7* (a, sin a) = j + . 
By (l2l?l) . h is decreasing in 7. So for any 7 > 7+ > 7*, we have 

h{a,M) < M«,/3,7 + ) < h(a,p,>y*) = 0. 

Combining this with (|26|) . we have To(tt,u) < h(a, /3, 7) < as desired. □ 

Proposition 7. 7/u is a pivotal point and 7 < 7 + , £/ien To(«,!i) < 0. 

Proof. We perform the following transformation. Move o n toward o„_i (the 
center of O n _i), along the x-axis while keeping et„_i and 6„_i on the boundary 
of O n until either 7 > 7+ or o n reaches o n _i. In the transformation, u stays 
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pivotal and the radius of O n changes. See Figure 2J The dotted circles show 
the process of the transformation. At the end of the transformation, if 7 > 7 + , 
then To(u,v) < by Proposition [51 and if o„ reaches o n _i, then Tq(u,v) = 
Tqj „_ 1 (u,f) < by the inductive hypothesis. Therefore, in order to prove 
that To(u,v) < before the transformation, it suffices to show that To(u,d) 
does not decrease during the transformation when X 0n , the x-coordinate of o„, 
decreases (i.e., we only need to show that dr °^ u ^ < q). 

Refer to Figure 0J Let X U ,Y U and X V ,Y V be the x- and y-coordinates of u 
and v, respectively. Let Y an _ 1 be the y-coordinatc of a n -i- Note that X U ,Y U , 
and Y a , n _ 1 remain constant during the transformation. Let r\ — f3r n . Since v 
stays pivotal during the transformation, r\ = (A n — B n )/2 = (\Po(u, &„-i)| — 
\Po(u, a n -i)\)/2 remains constant during the transformation. 

In the following, we express all other parameters as functions of (X Qn , 77, X Ul 
Yu, Y an l ) and then calculate their partial derivatives with respect to X Qn . This 
is possible because as mentioned above, 77, X u , Y u , and Y an _ 1 are independent 

y X 

of X n . Note that sin a = — — - and cos a = — . 



dr n d^Xl+Y^ Xon X 

— — — — — cos a. (33) 



da 9 W 2 + arctan(^-)) y y 



O-ri — 1 a n — l 



dx 0n dx 0n xi +y a 2 



(34) 



d(ar n ) dr n da 

— - — =a-^-- h r n — — = sin a - a cos a. (35) 

oX 0ri oX 0n oX 0n 

d(3 d(rj/r n ) r] dr n /3cosa 



dX Qn dX Dn r\ dX 0n r n 

dX v d(X 0n +r n cos(3) 80 dr n 
— =1 — r n smp— — hcosp- 



(36) 



dX 0n 0X On dX 0n dX 0n 

= 1-/3 cos a sin (3 — cos a cos f3. (37) 

dY v d(-r n sin f3) d/3 dr n 
- -r n cos Pttt^ sin p- 



dX 0n dX Qn dX 0n dX Qn 

= — /3 cos a cos /3 + cosasin/3. (38) 

Because Pq {u, and r\ remain constant during the transformation, from 
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([55)) we have 

d\Po(u,v)\ = d{P {u 1 a n - 1 ) + \A n \) 
dX 0n 0X On 
_ d\A n \ 
dX Dn 
_ d{ar n + rj) 

dX Dn 
_ d{ar n ) 
dX Qn 

= sin a — a cos a. (39) 

From Fact [5J H n — ||o n o„_i|| = X 0n — X 0n l . Since X 0n l is constant 
during the transformation, 

dH n = dXo^ = 
dX 0n dX Qn 

The vertical distance traveled by is r n — Y anl and the vertical distance 
traveled by Q^Li is r n ^\ — Y anl . When < a < ir/2, is inside O n _i and 
hence is red. So V n = (r„_i - Y an _ x ) - (r„ - l^-i) = r„_i - r n . Since r„_i is 
constant during the transformation, 

dV n dr n 

dxZ = -dxZ = cosa = lcosal 

When it/2 < a < n, both and Q^ n L 1 are green. So V n — (r„_i — Y an _ 1 ) + 
(r n — Y an l ). Since r„_i and Y an _ 1 are constant during the transformation, 

dV n dr n 

dx7rdxzr~ cosa= ^ cosaV 

Therefore, we have 

dSp _ dj^n ~ ri) - f Er=2 1 (2g, + ^)) 

_ip dX 0n 3 l «9x 0n ax J 

= — (/?cosa — |cosa|. (40) 

o o 
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Note that sin 7 = nfc^ and cos 7 = nf^jn^ • We have 

d\Do(u,v)\ = d\\uv\\ = d^/{X v - X u f + (Y v - Y u y 
dX Qn 0X On dX Qn 

_ X v — X u d(X v — X u ) Y v — Y u d(Y v — Y u ) 



\\uv\\ dX Qn \\ uv \\ dX 0n 

dX v . dY v 

cosl dxZ- smi dxZ (41) 

cos 7(1 — /3 cos a sin (3 — cos a cos j3) 
— sin 7(— (3 cos a cos (3 + cos a sin 0) 
cos 7 — cos a(/3 sin (3 cos 7 — j3 cos /3 sin 7 
+ cos /3 cos 7 + sin j3 sin 7) 

■ cos 7 — cosa(cos(/3 — 7) + /3sin(/3 — 7)). (42) 



Define a function 



fC fl ^ ,d\D {u,v)\ 



A(cos7 — cosa(cos(/? — 7) + /3sin(/3 — 7))). (43) 



From (El)— E3l), we have 



gTg(?V») = d\P (u,v)\ _ d\D (u 7 v)\ d\$o\ 

dX 0n dX Qn 0X On dX 0n 

lip tp 
= sin a — a cos a <p cos a — — | cosa| + /(a, (3, 7) (44) 

In the following, we will bound 9T d ° x ^' v} by single-variant functions in a. 

Refer to Figure SJ Let v" be the location of v when X Qn increases by dX 0n 
(i.e., o n moves to the right by dX 0n ). Let dl be the distance between v and v" . 
Let w be the angle from x-axis to vv . So cos a; = and sinw = t^-, where 
X v and 1^, are the x- and y-coordinates of v. Recall that 7 is the angle from uv 
to the a;-axis. Let 9 = u + 7. 
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From (|4ip . we have 



d\D (u,v)\ dX v dY v 



dX v dY„, dl 



cos 



dl ' dl ' dX 0n 

(cos 7 cos u> — sin 7 sin ui) 

cos(w + 7) 



dX 0n 

dl 



dX 0n 

= cose wz- ^ 

Therefore, / = -A a| ^fr' t ' )l = -A cos 9-^. Because gjt- > is independent 
of 7, when a and /3 are fixed the maximum of / occurs when cos# is minimized, 
i.e., when 9 is minimized or maximized. Since 9 = w + 7 where uj is independent 
of 7, 9 is minimized when 7 = and is maximized when 7 = 7+. Therefore 

/(a, /3, 7) < max{/(a, 0, 0),f(a, 0, 7 +)}, (46) 

where 7+ = 3sin , a ~" + arcsin ( - ?+ a j;° ) . 

4 V4Asin(— ^ )J 

On the other hand, || = A/3 cos a cos(/3 — 7). Since i> is an exit point of uv 
on the boundary of O n , —ir/2 < f3 — 7 < tt/2 and hence cos(/3 — 7) > 0. We 
distinguish two cases: 

1. tt/2 < a < tt. In this case, cos a < and as mentioned above cos(/3 — 7) > 
0. So || > when /3 < and || < when /3 > 0. This means that 

/(a,/3, 7 )</(a,0, 7 ). (47) 

When 7r/2 < a < it, we have |cosa| = — cos a, which, combined with 
(144|) . (|46|) . and gT]), implies that 



3To(w, v) 2ip 2ip 

— — — < sin ol — ol cos a — cos a 

aX 0n 3 3 

+ max{/(«,0,0),/(a,0, 7 + )}. (48) 



2. < a < tt/2. In this case, cos a > and as mentioned above cos(/3 
|| < when < 0, and % 



7) > 0. So || < when (3 < 0, and || > when /3 > 0. Observe 
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that /(a, /3, 7) — /(a,—/?, 7) = 2A cos a sin 7 (sin /3 — /3cos/3). By (|2"2"j) , 
< 7 < 7r/2 and hence sin7 > 0. Also sin0 - ft cos/3 > for < j3 < 1. 
Thus /(a, /3, 7) > /(a, -/3, 7) for < f3 < 1. This means that 

/(a,/?, 7) < /(a, |/3|, 7 ). (49) 

Since |i > in the range [|/3|,sina] and from (|2T|) . < sina, we have 

/(a, |0|,7) < /(a, sina, 7). (50) 

From (g!]) and (|50]) . 

/(a, 0, 7) < /(a, sin a, 7). (51) 



When < a < 7r/2, we have | cosa| = cos a, which, combined with (|44j) , 
(|4*5|) . andl5Tl implies that 

dT (u,v) 2tp Aif 

— < sin a — a cos a cos a 

dX 0n ~ 3 3 

+ max{/(a, sin a, 0), /(a, sin a, 7 + )}. (52) 

By (|4"5)) and (f52"|) , we only need to verify the following four inequalities. 

2tp 2tp 

31(a) = sin a — a cos a — cos a + /(a, 0, 0) 

o o 

< 0, when n/2 < a < it. (53) 
92(a) = sina — a cosa ^ ^ cosa + /(a, 0, 7 + ) < 

< 0, when ir/2 < a < n. (54) 
33(a) = sina — a cosa ^ ^ cosa + /(a, sina, 0) 

< 0, when < a < tt/2. (55) 

^ 2^ 4v? , ,/ • +\ 

.94(a) = sma — a cos a — cosa + / (a, sin a, 7 ) 

< 0, when < a < tt/2. (56) 

Since 31,32,53 and 34 are smooth functions on small intervals of a, one can 
easily verify the above inequalities using a numerical computing software, such 
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as Mathematica. We give a more formal approach in the appendix. We show 
that <?i,<72)53 and 54 have small Lipschitz constants, and then use a program 
that implements a simplified Piyavskii's algorithm |16j to verify that their upper 
bounds are less than 0. □ 

This completes the proof of Lemma [TJ □ 



5 Proof of Lemma 2 

This section contains the proof of Lemma [5J 

(proof of LemmaWi)- Let O be a set of chains whose stretch factor is greater 
than or equal to a threshold r. Suppose O is not empty. Let E be the subset 
of O consisting of chains with minimum number of circles; let it be n. For any 
chain O S E, associate with it a pair of terminals u 1 v with the largest stretch 
factor among all pairs of terminals of O. Fix a coordinate system. Without 
loss of generality, assume that all chains in E are normalized such that u is 
at the origin (0,0) and v is at (1,0). Each O £ E is represented by a vector 
x € lZ 3n that specifies, for each of the n circles in O, its radius and the x- 
and y-coordinates of its center. Since the conditions of the chain (Definition [1} 
and the conditions of E all include the boundary cases, E is represented by a 
non-empty closed set in lZ 3n . Define a function H (x) = Yyi=i r i- Since H(x) is 
continuous and has limit oo for ||x|| — > oo, H has a global minimum (see [111 
p.60]). Let O* be the chain in E that achieves the minimum of H. Then O* 
satisfies three conditions: (1) O* has stretch factor > r, (2) among all chains 
in O, the number of circles in O* is minimized, and (3) among all chains in E, 
the sum of the radii X^O so* r * ^ s minimized. 

Let u, v be the terminals associate with O* that yield the worst stretch 
factor. The chain O* has the following properties. 

Proposition 8. u and v are unobstructed in O* . 
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Proof. Suppose that Dq*(u,v) contains a point pj which is either cij or bj, 
1 < j < Ti — 1. Consider two sub-chains of O*: 0*j = (Oi, . . . ,Oj) and 
0* +ln — (Oj+i, ■ • ■ , O n ). So u,Pj are terminals of 0*j and Pj,v are ter- 
minals of 0* +l n . Note that \Po*(u,v)\ < \Pq* . (u,pj)\ + \Po* +ln (Pj,v)\ and 
\D *(u,v)\ = \D * .(u,pj)\ + \Do» +ln (pj,v)\. So we have either \P * .(u,Pj)\ > 
p\Do* .(u,pj)\ or |-Pe>? +1 (Pj, v )\ > p\Do* j (pj,v)\. This means that either 
0\ ■ or 0* +1 n has stretch factor greater than or equal to that of O* and has 
less number of circles than O* — a contradiction to condition (2) of O* . So u,v 
must be unobstructed in O* . □ 

Proposition 9. B1B2 ■ ■ ■ B n is a shortest path between u and v in O* and so 
is A1A2 . . . A„ . 

Proof. Suppose that P>\P>2 ■ ■ ■ B n is not a shortest path between u and v. Let Bj 
be the first arc in the path B\Bi ■ ■ ■ B n that is not in a shortest path between u 
and v. This means that every shortest path between u and v contains Aj. This 
also means that Bj-i and aj—ibj-i belong to a shortest path between u and 
v. We can assume that Bj is not degenerated because otherwise ajbj will be a 
shortcut between bj-\ and aj in the shortest path containing Bj-%, aj—ibj—i 
and Aj. 

We perform a transformation that shrinks Oj as follows: keep Oj_i and aj 
on the boundary of Oj and reduce the radius Tj by a small amount; in case 
where Aj is degenerated, keep aj {aj~i) on the boundary of Oj and move the 
center Oj by a small amount along the line segment ojaj toward aj. See Figure [5] 
for an illustration. Such transformation can be performed while maintaining the 
following: 

1. O* remains a chain. Since aj-i and aj stay on the boundary of Oj, Oj 
remains intersected with Oj-i and 0_,+i during the transformation. So 
property (1) of chain (in Definition [1]) is satisfied. Shrinking Oj will shrink 

(?) (i) (i—2) (7+2) 

the connecting arcs Cyj_ x and C)Xi, but does not affect C) J _ X ' and C) J + X '. 

(7—2) (?} (j) 

So Cj_ x and Cj_ x do not overlap after shrinking Oj, and neither do C^+i 
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Figure 5: Illustrations of Propositions |H1 shrinking a circle. Two cases are 
shown. Aj is not degenerated in case 1 and Aj is degenerated in case 2. 

and Cj^ 2 '. The only other possible way that the property (2) of chain 
can be violated is for C^ 1 ' and C^ +i ' to overlap; but we know that Bj 
is not degenerated. So shrinking Oj by a sufficiently small amount will 
not make C^~^ and C^ +r> overlap. 

2. Bj is still not in a shortest path between u and v. This is because the 
lengths of the paths between u and v changes continuously when Oj 
shrinks. So we can shrink Oj by a sufficiently small amount such that 
Bj is still not in a shortest path between u and v. 

3. \Dq*(u,v)\ remains unchanged. This is because by Proposition [51 u,v are 
unobstructed in O* . So we can shrink Oj by a sufficiently small amount 
such that u, v are still unobstructed, which means that \Dq*(u,v)\ = \\uv\\ 
remains the same. 

We claim that shrinking Oj by this transformation will not decrease \Po* (u, v) | . 
Here is why. Shrinking Oi does not affect O* except for Aj and the arcs and 
line segments containing bj—\ or bj, including -Bj-i, Bj, Bj + \, a,j-\bj-i, and 
ajbj. 

First consider the effect of shrinking Oj on \Aj\. Since Bj is not in a shortest 
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path, \Aj\ < ||a_y_i6j_i|| + \Bj\ + \ \a,jbj\\ and hence \Aj\ is less than half of the 
perimeter of Oj. Therefore, shrinking Oj by a sufficiently small amount while 
keeping Oj-i and aj on the boundary of Oj will increase \ Aj\. 

Next consider the effect of shrinking Oj on a shortest path containing bj-\ 
or bj. Let P be a shortest path between u and v in O* containing bj. Since 
P does not contain Bj, P includes Aj, ajbj, and -Bj+i- Let b'j be the new 
location of bj and let bjb'j be arc between bj and 6^ that becomes an extension 
of -Bj+i after the transformation. So shrinking increases |-Bj+i| by \bjbj \ and 
decreases 1 1 aj bj\\ by 1 1 aj bj\ \ — \ \ aj b'J\. Observe that 

l¥?l > IMjH - llo^H- (57) 

For example, in Figure [5j b' 4 is the new location of 64 after shrinking O4 and 
I6464I > 1 1 «4 £»4 1 1 — I ) 1 1 ■ This implies that shrinking Oi increases |-Bj+i| + 
||aj6j||; i.e., shrinking Oi increases the length of P. By a similar argument, if 
P contains bj-\, shrinking Oi increases the length of P. 

So in any case, shrinking Oj will not decrease \Pq* (u, v)\. 

In summary, by shrinking Oj we have a new chain that satisfies condition (1) 
and (2) of O* , and has a smaller sum of radii than O* , which is a contradiction 
to condition (3) of O* . 

This proves that B1B2 ■ ■ ■ B n is a shortest path between u and v in O* . 
Similarly, A\Ai . . . A n is also a shortest path between u and v. This completes 
the proof. □ 

Now we are ready to prove that $0* > — ■^^■\Po* (u, v)\. 

See Figure [5] We align the centers 01 . . . o n of the circles in O* on a straight 
line I by the following operations: first rotate 0\ around o\ such that u comes 
to a location where \A\\ ~ \B\\; then rotate 0\ and O2 around 02 such that 
|^2| = I-B2I; then rotate Oi,02, and O3 around 03 such that \A S \ — \B 3 \; and 
so on. Let O* be the resulting chain after the rotation operations. In O*, we 
have A{ — B t for 1 < i < n and X^r=i(l^i| + 1^1) in C* remains the same as in 
O*, which is 2\Pq* (u, v)\ by Proposition |H1 So the length of the path Ai ...A n 
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Figure 6: The centers of the circles in O* are aligned on a straight line. The 
path V-qt in which the overlapping portions of red arcs and green arcs cancel 
each other. 

in O* is 

|^ 1 ...i4 n | = |P .(«,«)|. (58) 

Recall that $o» = (p(r n - n) — | ^™ =2 (2#j + Vi). By reversing the labels 
of the circles in O* , if necessary, we can assume that r n >r\. Now it suffices 
to show that Ya=2^ h i + Vi) < V5\Po* (u, v)\. 

Recall from Definition [5] that Hi and V% are the horizontal and vertical dis- 
tance traveled along the path Vi, where Vi is a path from q~L 11 a peak of Oj-i, 
to qf , a peak of Oi. Also recall that Vi consists of green or red arcs, the green 
arcs in Vi contribute positively to Hi and Vi, and the red arcs in Vi contribute 
negatively to Hi and Vi . Since the centers o\ . . . o n are all aligned on a straight 
line in O* , the two peaks and q~* of every circle Oi overlap, for 2 < i < n — 1. 
So the paths 7^ are joined at the peaks to form a single path from q^ to q^ ; 
denote it by V-qt. Refer to Figure El Let the overlapping portions of red arcs 
and green arcs in V-qt cancel each other out and then remove the remaining red 
arcs, if any, at the beginning or the end of V-^r. The resulting path, denoted 
by V^, is a subpath of A\ . . . A n . Thus from ([51]) 

\VL F \<\A l ...A n \ = \P «{u,v)\. (59) 

Let H- — and V— be the horizontal and vertical distance traveled by V- — . By 
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the definition of VL-, J2?=2 Hi < and £™ =2 V t < V^. Combining these 
with (|59"1) . we have 



(E + (E ^ V(^) 2 + (^) 2 

=2 »=2 

<|^r|<|ib.(u,t;)|. (60) 



It follows that 



n 

2 



2 E^ + E ^ ^ J 5( ^ + ^ 

i=2 i=2 y i=2 i=2 

<V5\P -(u,v)\, (61) 
as required. This completes the proof of Lemma 2. □ 



6 Conclusions 

In this paper, we showed that the stretch factor of the Delaunay triangulation 
is less than 1.998 by proving the same upper bound on the stretch factor of a 
chain. 

There are a few places where our approach can be further improved. First, 
the potential function can be improved to yield a better upper bound. For ex- 
ample, if we define the potential function <&q to be the length of the segment 
of uv inside O n , then we can improve the upper bound to 1.98, although the 
analysis is quite complicated. Secondly, the key components of our approach are 
the proofs of Proposition [6] and Proposition [Jj which rely largely on functional 
analysis. We hope to gain insight of the underlying geometry that will help us 
simplify the proofs and push the upper bound closer to the tight bound, which 
we believe is around 1.6. 
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7 Appendix 

7.1 Verify Inequalities fl55]> . (|54j), (|55|), and (|56l) 

Let 2 = a/4 + sina/4. Since < z < 7r/4, we have 1 < < 1.111. Further- 
more, 



dz d(a/4 + sina/4) 



da 



da 



= 1/4+ cos a/4, 



and 



1 2 COS z 1 cosz , „ ,„„ 

< : = tan(z/2) < 0.415. 



dz sinz (sinz) 2 — sinz sinz 
Recall that 

3 sin a — a 



7 
So 



a + sm a 



3 sin a — a . . z 

, — h arcsml ). 

4Asin(^i±f^) ) 4 v Asinz ; 



d( 7 H 



da 



< 



3 cos a — 1 



1 



d(-M dz 

V sin 2 / 



< 1 



< l 



4 A,/l- dz da 

V ^ Asms'' 

°- 415 

= (1/4 + cos a/4) 



Av/I - (t-?— V 

V v A sin 2 / 

0.415 



(1/4 + 1/4) 



< 1.15. 

Recall that 

f(a,P, 7) = — A(cos7 — cosa(cos(/? — 7) + /3sin(/3 — 7))). 
We have 

df(aJ)A)) d(-\h -msMl 

A sin a < 2. 



(62) 



df (a, 0,0) 




d(— A(l — cos a)) 




da 




da 
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rf/(«,0,7 + ) 
da 



d(— A(cos7+ — cosacos(— 7+))) 



da 

d(— A cos 7+ (1 — cos a)) 



< 1.8(1.15 + 1) < 4. 



sin7 + (l — cos a) | + | sin a cos 7+ 1) 



df(a, sin a, 0) 




A(l — cosacos(sina) — cosasinasin(sino;))) 


da 







< A(| sinacos(sina)| + |(cosa) 2 sin a cos (sin a) \ + |(sina:) 2 sin(sina)|) 

< 1.8(1 + 1 + 1) < 6. 



df(a, sin a, 7+) 



da 



d{— A(cos7+ — cos a cos (sin a — 7+) — cosasinasin(sina! — 7 + ))) 



Also 



da 

d{— A(cos7+ — cos a cos (sin a — 7+) — Sln ( 2a ) sin(sina: — 7 + ))) 



da 
d(-f + ) 

< A(| — - — sin7 + | + I sinacos(sina — 7 + )| 
da 

1 • , ■ -i-N/ d(-f + ).. 

+ I cos a sin(sin a — 7 T )(cosa — ) 

da 

...... ... sin(2a) , rf(7 + ) 

+ I cos(2a) sm(sma — 7 )| + | ^ cos(sma — 7 T )(cos a )|) 



da 



< 1.8(1.15 + 1 + (1 + 1.15) + 1 + (1 + 1.15)/2) < 12. 



d(sin a — 


a cos a — 


2ip 
3 


-j£ cos a) 






da 








d (sin a — 


a cos a — 


%£_ 

3 


^ cos a) 




da 





2w 2ip 
\a sin a + — sin a\ < tt + — < 4. 
o o 



Aw iip 
\a sin a + — sin a\ < tt + — < 4. 
3 3 
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Therefore the Lipschitz constants are 



dgi(a) 




d(sin a — 


a cos a — 


2t£ 

3 


2^ 
3 


cos a) 


da 






da 








dg 2 (a) 




e?(sin a — 


a cos a — 


2ip 
3 


3 


cos a) 


da 






da 








dg 3 (a) 




d(sin a ~ 


a cos a — 


2ip 
3 


3 


cos a) 


da 






da 








dg 4 (a) 




e?(sin a — 


a cos a — 


2ip 
3 


2^ 
3 


cos a) 


da 




da 



d/(a,0,0) 



+ 



tfa 

#(a,0, 7 + 



d/(a, sin a, 0) 



da 

df{a, sin a, 7 + ) 



< 4 + 2 < 16. 
< 4 + 4 < 16. 
< 4 + 6 < 16. 



da 



< 4 + 12 = 16 



Now we can use a simplified version of the Piyavskii's algorithm |16j for 
Lipschitz optimization to verify that gt(a) < 0, 1 < i < 4. The following 
Algorithm Bound((^, s, t) will either find a value gi(a) > in the given range 
or return an upper bound on the value of gi in range [s,t] that is less than 0. 
We run Bound(<^, s, t) with s and t set to be the lower and upper bound on 
the range of a and verify that the upper bound of gi(a) is indeed less than in 
that range, for 1 < i < 4, where gi(0) is set to be lim Q _j.Q gt{a). 



Algorithm Bound(g;, s,t), for 1 < i < 4 



1. if gi(s) > or gi(t) > return max{<?;(s), gi(t)}. 

2. apex = ma 1 x{gi(s) , gi(t)} + 16 * (t — s)/2 

3. if apex > 0, do: 

3.1. apex = max{Bound(</i, s, (s + t)/2), Boundfji, (s + t)/2, t)} 

4. return apex 
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